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We calculate the time evolution of a far-from-equilibrium initial state of a non-relativistic 
ultracold Bose gas in one spatial dimension. The non-perturbative approximation scheme 
is based on a systematic expansion of the two-particle irreducible effective action in powers 
of the inverse number of field components. This yields dynamic equations which contain 
direct scattering, memory and off-shell effects that are not captured in mean-field theory. 

1. Non-equilibrium ultracold atomic gases 

In recent years, sophisticated new technologies have been developed for the cooling 
and trapping of ultracold atomic gases: These allow the contact-free storage of atoms 
in almost arbitrarily shaped "atom traps", as well as the manipulation and Feshbach- 
resonant enhancement of the interactions between the atoms using externally applied 
electromagnetic fields [[Q. (Quasi) one- and two-dimensional traps [0E] as well as optical 
lattices [ @j allow to realize strongly correlated many-body states of atoms reminiscent 
of similar phenomena in condensed matter systems. Since these external conditions can 
be varied very quickly while the experiment is running and the reaction of the gas be 
monitored extremely precisely, high-precision studies of the quantum dynamics of many- 
body systems become possible even far from equilibrium. 

The 1+1-dimensional systems of identical bosons of mass m to be considered are de- 
scribed by fields obeying [^f(x, t), $?'(y, t)] = 5(x — y). Their dynamics is governed by the 
Hamiltonian 

H = - dxVUx) -tf(a;) + - dxdyVUx)VUy)V(\x-y\)V(y)V(x). (1) 

J 2m 2 J 

At ultralow temperatures, where the thermal de Broglie wave lengths are much larger than 
the effective range of the binary Born-Oppenheimer potential V, this can be approximated 
by the effective local coupling V{x) = gmS(x). 

A particular challenge represent systems far from equilibrium for which conventional 
perturbative and mean-field approximation schemes fail to describe the long-time evolu- 
tion. We present a dynamical many body theory of an ultracold Bose gas which system- 
atically extends beyond such approximations. The non-perturbative approach is based on 

* Funded by Deutsche Forschungsgemeinschaft. 



2 



T. Gasenzer, J. Berges, M. G. Schmidt, and M. Seco 



a systematic expansion of the two-particle 
irreducible (2PI) effective action in powers 
of the inverse number of field components 
M [0 0Hj. This 2PI 1/7V expansion, 
to next-to- leading order (NLO), yields 
dynamic equations which allow to describe 
far-from-equilibrium dynamics as well as 
the late-time approach to quantum thermal 
equilibrium. We emphasize that mean-field 
approximations fail to describe this dynam- 
ics even qualitatively Similarly, standard 
kinetic descriptions based on two-to-two 
collisions give a trivial (constant) dynamics 
because of phase space restrictions in one 
spatial dimension. Recently, these methods 
have allowed important progress in describ- 
ing the non-equilibrium dynamics of strongly interacting relativistic systems for bosonic 

[El El El HHj, as well as fermionic degrees of freedom [|TTJ [T2]. The 2PI I/TV expansion 
has also been successfully applied to compute critical exponents in thermal equilibrium 
near the second-order phase transition of a model in the same universality class [115]. For 
details we refer to [ Hlj as well as to [ IT5| ITo]. 

2. The 2PI effective action in NLO l/Af approximation 

The 2PI effective action [ El IS] is obtained as a double Legendre transform of the 
generating functional of connected Greens functions and is conveniently written as 

T[<j>, G] = S[(f)] + % - Tr {lnC- 1 + G l [(f)]G} + r 2 [0, G] + const. , (2) 

which contains the contribution from the classical action S, a one-loop-type term and a 
term T 2 [(f), G] that contains all the rest. Gq 1 is the inverse of the classical propagator with 
iGq y\ (p) = 8 2 S[(f>}/8(j>i(x)8<j>j{i)). We denote the time-space vector by x = (x ,Xi), 
etc., where Xq is the time component. 0,, i = 1, 2 are the Af = 2 components of the macro- 
scopic or mean field, e.g., the real and imaginary parts (j>% = y/2Ke 2 = ^/2Im Cor- 
respondingly, the components of the connected 2-point function are denoted by Gij(x, y) = 
l(T<f)i(x)<f)j{y)) — (pi(x)(pj(y). The equations of motion for <pi{x) and Gij(x,y) are, finally, 
derived using the stationarity conditions ST[4>, G]/8(j)i{x) = 0, ST[(f>, G]/5Gij(x, y) = 0. 
T 2 may be represented by the series of all closed 2PI diagrams consisting of full prop- 
agators G, field insertions 0, and bare vertices proportional to #id [ 5j. The 2PI 1/jV 
expansion involves a resummation of the diagrams up to infinite order in gxo and reaches 
substantially beyond the conventional (1PI) l/Af approach. For the NLO resummation 
cf . Fig. ^ for explicit expressions of T and the dynamic equations we refer to [ [Hi • For 
the following discussion, it is convenient to express all 2-point correlation functions in 
terms of their connected statistical and spectral parts, Fij(x,y) = $j(y)}) c /2, 
pij(x,y) = i([$i(x),$j(y)]) c , such that Gij(x,y) = Fij(x,y) - (i/2)p ij (x, y) sgn c (x - y ), 
where C denotes ordering along the closed real-time path. 




Figure 1. Diagrammar of the LO and NLO 
contributions in the l/AA-expansion, to the 
2P1 part T 2 [(j), G] of the 2PI effective action. 
The solid lines represent 2-point functions 
Gij(x,y), the crosses fields (fii(x), and the 
wiggly lines vertices gir>8(x — y). 
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Figure 2. Equilibration of ultracold ID Bose gas. (a) and (b) Time evolution of momen- 
tum mode distribution, (c) Development of thermal distribution with temperature G. 
(inset) Momentum dependence of self energy, (d) Establishment of fluctuation-dissipation 
relation, (inset) Energies. For details confer Section 

3. Equilibration of a homogeneous ultracold 1-dimensional (ID) Bose gas 

We have investigated the dynamic evolution of a ID Bose gas of sodium atoms in a box 
of length L = N s a s (N s = 64 modes on a grid with a s = 1.33 /zm), with periodic boundary 
conditions, starting in a non-equilibrium state which is solely described by a Gaussian 
momentum distribution n(t — 0;p) — N~ l exp{— p 2 /a 2 } and normalized such that the 
line density is fixed to n\ = 10 7 atoms/m. All higher correlation functions as well as <pi 
are assumed to vanish initially [E]- The atoms are weakly interacting with each other, 



such that #id = h 2 'jni/m, with the dimensionless parameter 7 = 1.5 ■ 10~ 3 . Figs. |2K ; b 
show the time evolution of the occupations of the momentum modes Pi = (2/a s ) sin(7ri/N s ) 
as functions of p and t, respectively. Although the system very quickly, after about 5 /zs, 
evolves to a quasistationary state, the final drift to the equilibrium distribution takes 
roughly ten times longer. Note that the mean-field Hartree-Fock (HF) approximation, 
which only takes into account T^ and the 0(giD)-diagram in Fig. ^p, would conserve 
exactly all mode occupations and no equilibration would be seen. We further investigated 
the different time scales by fitting the distribution to the Bose-Einstein form n(t;p) = 
[exp{(p 2 /2m — /i)//cs0(t;p)} — with a p-dependent temperature variable Q(t;p). 
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Fig. Efc shows 0(t;p) for t = 0...0.6s. Hence, during the quasistationary drift period, no 
temperature can be attributed to n(t;p), while, for large t, © becomes approximately p- 
independent. Its variation at small p reflects the modified dispersion relation. At large p, 
excitations of the gas will be approximately single-particle-like, such that we can associate 
0(0.6 s; 128/L) with the final temperature T. Using this, we deduce the shift in the 
dispersion relation, n(0.6s;p) = [exp{(p 2 /2m + — fi)/k B T} — (inset of Fig. |2t) 
and, with the total energy in equilibrium, E tot = ^2 \p 2 /2m + E(p)/2]n(0.6s;p), which 
must be equal to the initial one, the chemical potential /i = 1.08 gm^i- 

We studied the time dependence of the ratio of unequal-time correlation functions, i.e., 
= [(F 11 (t,0;p) 2 + F 12 (t,0;p) 2 )/( Pll (t,0;p) 2 + p 12 (t,0;p) 2 )] 1 /Vn(t;p), Fig. EH. £ is 
a measure of the interdependence of the statistical and spectral functions, which, in ther- 
mal equilibrium, are connected through the fluctuation-dissipation relation F^ eq \to,p) = 
—i[n(u, T) + l/2]p( eq )(uj,p). One finds that already during the quasistationary drift period, 
£ becomes approximately time-independent for all momenta Pi, i.e., there is a fluctuation- 
dissipation relation even though the system is still away from equilibrium. A similar 
signature is found when comparing the kinetic and interaction contributions to the total 
energy as shown in the inset of Fig. |2Ji. During the drift period, these contributions are 
constant and approximately equal to each other, calling in mind the virial theorem. In 
summary, during the drift period, the system is not yet in equilibrium as far as the mo- 
mentum distribution and temperature is concerned, but shows important characteristics 
of a system close to equilibrium. The details of this is a subject of future work. 
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